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Abstract — We study the asymptotic properties of distributed 
consensus algorithms over switching directed random networks. 
More specifically, we focus on consensus algorithms over inde- 
pendent and identically distributed, directed random graphs, 
where each agent can communicate with any other agent with 
some exogenously specified probability. While different aspects 
of consensus algorithms over random switching networks have 
been widely studied, a complete characterization of the distri- 
bution of the asymptotic value for general asymmetric random 
consensus algorithms remains an open problem. In this paper, 
we derive closed-form expressions for the mean and an upper 
bound for the variance of the asymptotic consensus value, when 
the underlying network evolves according to an i.i.d. directed 
random graph process. We also provide numerical simulations 
that illustrate our results. 

I. Introduction 

Distributed consensus algorithms have attracted a signif- 
icant amount of attention in the past few years. Besides 
their wide range of applications in distributed and parallel 
computation [1], distributed control [2], [3] and robotics [4], 
they have also been used as models of opinion dynamics and 
belief formation in social networks [5], [6]. The central focus 
in this vast body of literature is to study whether a group of 
agents in a network, with local communication capabilities 
can reach a global agreement, using simple, deterministic 
information exchange protocols^ 

In recent years, there has also been some interest in un- 
derstanding the behavior of consensus algorithms in random 
settings [7]-[12]. The randomness can be either due to the 
choice of a randomized network communication protocol 
or simply caused by the potential unpredictability of the 
environment in which the distributed consensus algorithm 
is implemented [13]. It is recently shown that consensus 
algorithms over i.i.d. random networks lead to a global agree- 
ment on a possibly random value, as long as the network 
is connected in expectation [10]. While different aspects of 
consensus algorithms over random switching networks, such 
as conditions for convergence [7]-[10] and the speed of con- 
vergence [13], have been widely studied, a characterization 
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of the distribution of the asymptotic consensus value has 
attracted little attention. Two notable exceptions are Boyd et 
al. [14], who study the asymptotic behavior of the random 
consensus value in the special case of symmetric networks, 
and Tahbaz-Salehi and Jadbabaie [15], who compute the 
mean and variance of the consensus value for general i.i.d. 
graph processes. Nevertheless, a complete characterization of 
the distribution of the asymptotic value for general asymmet- 
ric random consensus algorithms remains an open problem. 

In this paper, we study asymptotic properties of consensus 
algorithms over a general class of switching, directed random 
graphs. More specifically, building on the results of [15], 
we derive closed-form expressions for the mean and an 
upper-bound for the variance of the asymptotic consensus 
value, when the underlying network evolves according to 
an i.i.d. directed random graph process. In our model, at 
each time period, a directed communication link is estab- 
lished between two agents with some exogenously specified 
probability. Due to the potential asymmetry in pairwise 
communications between different agents, the asymptotic 
value of consensus is not guaranteed to be the average of the 
initial conditions. Instead, agents will asymptotically agree 
on some random value in the convex hull of the initial 
conditions. Furthermore, our closed-form characterization 
of the variance provides a quantitative measure of how 
dispersed the random agreement point is around the average 
of the initial conditions in terms of the fundamentals of the 
model, namely, the structure of the network, the exogenous 
probabilities of communication, and the initial conditions. 

The rest of the paper is organized as follows. In the 
next section, we describe our model of random consensus 
algorithms. In Sections [HI] and [IV] we derive an explicit 
expression for the mean and an upper bound for the variance 
of the limiting consensus value over switching directed 
random graphs, respectively. Section V contains simulations 
of our results and Section VI concludes the paper. 

II. Consensus Over Switching Directed Random 
Graphs 

Consider the discrete-time linear dynamical system 
x(k)=W k x(k-l), 



(1) 



where lc G {1,2,...} is the discrete time index, x(k) 6 M" 
is the state vector at time lc, and {Wt}~ =1 is a sequence of 
stochastic matrices. We interpret (HJ as a distributed scheme 
where a collection of agents, labeled 1 through n, update 
their state values as a convex combination of the state values 
of their neighbors at the previous time step. Given this 
interpretation, x,-(fc) corresponds to the state value of agent i 



at time k, and Wk captures the neighborhood relation between 
different agents at time k: the ij element of Wk is positive 
only if agent i has access to the state of agent j. For the 
remainder of the paper, we assume that the weight matrices 
Wk are randomly generated by an independent and identically 
distributed matrix process. 

We say that the dynamical system (Q]) reaches consensus 
asymptotically on some path {WIt}™ = j, if along that path, 
there exists x* G R such that Xj(k) —> x* for all ; as k — > °°. 
We refer to x* as the asymptotic consensus value. It is well- 
known that for i.i.d. random networks, the dynamical system 
(Q]l reaches consensus on almost all paths if and only if the 
graph corresponding to the communications between agents 
is connected in expectation. More precisely, it was shown 
in [10] that Wk---WiWi — > ld T almost surely — where d 
is some random vector — if and only if the second largest 
eigenvalue modulus of EWk is subunit. Clearly, under such 
conditions, the dynamical system in (Q]) reaches consensus 
with probability one where the consensus value is a random 
variable equal to x* = d T x(0), where x(0) is the vector of 
initial conditions. 

A complete characterization of the random consensus 
value x* is an open problem. However, it is possible to 
compute its mean and variance in terms of the first two 
moments of the i.i.d. weight matrix process. In [15], the 
authors prove that the conditional mean of the random 
consensus value is given by the random consensus value are 
given by 



Ex* = x(0yv l (EWk), 



(2) 



and its conditional variance is equal to 

var(x*) = [x(0)®x(0)fvec(cov(d)) (3) 
= [x(0)®x(0)] T Vi(E[W*®W*]) - [x(0) T \i(EW k )] 2 

where Vi (■) denotes the normalized left eigenvector corre- 
sponding to the unit eigenvalue, and (g> denotes the Kronecker 
product. In [16], the authors derive an explicit expression 
for the mean and variance of x* in the particular case 
of a switching Erdos-Renyi random graph process. In the 
following, we shall use ([3]) to extend these results and derive 
an explicit expression for the mean and an upper-bound for 
the variance of the asymptotic consensus value over a wider 
class of switching, directed random graph processes. 

III. Mean Analysis for Directed Random Graphs 
Processes 

A. Directed Random Graph Process 

Consider a connected undirected graph G c = {V,E C ) with 
a fixed set of vertices V = [«] , and unweighted edges (no 
self-loops allowed). Each undirected edge in E c represents 
a potential communication channel between nodes i and j, 
where this channel can be used to send information in both 
directions. In this paper, we focus on directed communica- 
tions, i.e., the event of node ; sending information towards 
node j is independent from the even of node j sending 
information towards i. In this context, it is convenient to 
interpret an undirected edge {i,j} G E c as the union of two 



independent directed edges, , (j, i)}, where the ordered 

pair represents a directed link from node i to node j. 

In this paper, we study randomized time-switching con- 
sensus processes. In particular, in each discrete time slot 
k > 1, we construct a random directed graph = (V, $k), 
with $k C E c , such that the existence of a directed edge 
(w,v) G §k is determined randomly and independently of all 
other directed edges (including the reciprocal edge (v,u)) 
with a probability p uv G (0, 1) for (u, v) G E c , and p uv = for 
(w,v) $ E p . In other words, in each time slot, we randomly 
select a subset §k of directed links chosen from a set of 
candidate (directed) links in E c . We are specially interested 
in the case in which the probability p uv of existence of a 
directed link (u,v) depends exclusively on the node that 
receives information via that link, i.e., Pr((«,v) G $k) = p v , 
where p Y G (0, 1). In this setting, we can model the ability of 
a node to 'listen' to their neighboring nodes. For example, 
in the context of opinion dynamics in social networks [5], 
[6], [12], the probability p v can represent the tendency of the 
individual at node v to take into account the opinion of her 
neighbors (which could depend, for example, on how many 
acquaintances the individual has). 

Let us denote by A c the symmetric adjacency matrix 
of the graph G £ , where entries ajj = 1 if {i.J} G E c , and 
a\j = otherwise. We define the degree of node i as di = 
YJ]=\ a jii and the associated degree matrix as D c = diag(di). 
We also denote the random (nonsymmetric) adjacency matrix 



associated with as ^k = 
as 



Jk) 



which can be described 



a (k) 



w.p. p v 
w.p. 1 - 



(4) 



We denote the in-degrees of as cffl = and the in- 

degree matrix as £>k = diag(dy). From the definition of Sf^, 
the in-degrees are independent Bernoulli random variables 
dy ~Ber(rf v ,p v ), i.e., 

We describe the consensus dynamics in ([T]i via a sequence 
of stochastic matrices {Wit}~ =1 associated to the sequence of 
random directed graphs {5^ = (y,$k)Yk=\ as follows: 

W k ={£>k + I)~ i {Ak + I)\ for*>l. (5) 

Notice that adding the identity matrix to the adjacency in 
(|5]l is equivalent to adding a self-loop to every vertex in V. 
We include these self-loops to avoid singularities associated 
with the existence of isolated nodes in SJJfc (f° r which di = 0, 
and £>k would not invertible). Since G c is connected and 
p v > for all v G V, the expected communications graph is 
connected and the stochastic dynamical system in (Q]) reaches 
consensus on almost all paths, although the asymptotic 
consensus value x* is a random variable (not the initial 
average). 

In the following subsections, we first derive closed-form 
expressions for the mean of x*, and an upper-bound for the 
variance var{x*) in terms of the following elements: 



1) The set of initial conditions, {x u (0)} ueV , 

2) the set of nodes properties, {(pu,d u )} ueV , and 

3) the network topology, via the eigenvalues of the ex- 
pected matrix EW k . 

As we shall show in Section |IV] our expression for the 
variance has a nice interpretation, since it separates the 
influence of each one of the above elements into three 
multiplicative terms. In the next subsection we provide the 
details regarding our analysis of the expectation of x*. 

B. Mean of Consensus Value 

We use (fJJ to study the mean of the consensus value. 
We first derive an expression for EW k , and then study its 
dominant left eigenvector vi (EW k ). For notational conve- 
nience, we define the random variable Zi = 1/ (di + l) where 
di ~Ber(di,pi), and denote its first and second moments as 
M ( (1) = E (zi) and = E (zj) . The diagonal entries of EW k 
are then given by Ewu = E [l/ (d{ + l)l , which present the 
following explicit expression (see Appendix I for details): 



Ewu=MP= 



(6) 



Pi (di + l) 

where = 1 — pi. Furthermore, the off-diagonal entries of 
EWk are equal to (see Appendix I for details): 

di+l 



Ewi 



1i 



+ Pi(di+l)-\ 



Pi(di+i)di 



1 -M; 

di 



(7) 



Taking (O and (0 into account, we can write EW k as follows: 

EW k =Z+(I-Z)D- 1 Al, (8) 



where E = diag 



M, 



i) 



As expected, it is easy to check 
that EWk is a stochastic matrix, i.e., (EW^) 1„ = El„ + 
{I-JL)D- l A T c l n = l n . 

Based on ([8]l we can write E(x*) explicitly in terms of di 
and pi, as follows: 

Theorem 1: Consider the random adjacency matrix A k in 
(O and the associated (random) stochastic matrix W k in (0. 
The expectation of the asymptotic consensus value of (Q~|) is 
given by 

E(x*) = £pw ( - */(()), (9) 
i=i 

where 

Pi(di+ 1) di 



Wi(pi,di) 



P(Pi4i 



Pi (di + l)-l 

Y< w i(Pi'di) 



di+l 



and 



Proof: Our proof is based on computing vi (EWj) 
and applying Let us define v = \'i (EW)t) and w = 
(/ — E) v. From ([8j, we have that the eigenvalue equation 
corresponding to the dominant left eigenvector of EW k is 



v T (E+(/-E)D- 1 A,^ 



which can be rewritten as 



\ T (I — L)D~ 1 AT = \ T (I — E). This last equation can be 
written as w T D^ l A^ =w T . The solution to this equation is 
the stationary distribution of the Markov chain with transition 
matrix D~ A^ , which is equal to n = d/^,rf,, where d = 
(di, ...,d„) T . Hence, the solution to the eigenvalue equation 
is 

V! (EW k ) = a{I-Y.y l d, 
where we include the normalizing parameter 

/ v- 1 di \ 



(l) 



\ / 1 -M\ 

such that (EWi)||j = 1. Hence, from (O we have 



~ 1-M, 



ci(0). 



Substituting the expression for in ©, we reach (0 via 
simple algebraic simplifications. ■ 

In general, the asymptotic mean E(x*) does not coincide 
with the initial average jco = -L; x /(0)- There is a simple 
technique, based on Theorem [T] that allows us to make the 
expected consensus value to be equal to the initial average. 
This technique consists of using y, (0) = (pwi) l Xi (0) as 
initial conditions in ([T). Hence, one can easily check that the 
asymptotic consensus value E(y*) equals the initial average 

X{). 



IV. Variance of the Asymptotic Consensus Value 

In this section, we derive an expression that explicitly 
relates the variance var(x*) with the three elements that 
influences it, namely, the set of initial conditions {x u (0)} ueV , 
the nodes properties {(p,i,d u )} ueV , and the network structure 
(via the eigenvalues of the expected matrix EWj;). For sim- 
plicity in notation, we denote E{W k ®W k ) and iEW k ®EW k ) 
by R and Q, respectively. Our analysis starts in expression 
(0, which can be rewritten as 

var (x*) = [x(0) ® x(0)] r [vi (R) - vj (Q)} . 

Hence, we can upper-bound the variance of the asymptotic 
consensus value as follows: 



•(af*)<||x(0)(g)x(0)|| 1 ||v 1 (R)-n(Q)\l 



(10) 



From the rules of Kronecker multiplication, we can write the 
first factor in terms of the initial conditions as 



|x(0)<g>x(0)| 



I Vi (0)^(0)| 



(11) 



1 <i,j<n 



In the following, we derive an upper bound for the sec- 
ond factor ||vi (/?) — vi ((2)|L in terms of the nodes prop- 
erties and the network structure. Our approach to bound 
|| vi (R) — vi (Q)!!^ is based on the observation that both R 
and Q are n 2 x n 2 stochastic matrices, and the dominant left 
eigenvectors Vi (R) and Vi (Q) are stationary distributions of 
the Markov chains with transition matrices R and Q. We 
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Fig. 1. The pattern of E[Wi53 Wj] for n = 3. The numbers in parentheses 
represent the value of each entry defined in the Appendices. 



denote these distributions by Vi (R) = % and vi (Q) = %, re- 
spectively. In this setting, we can apply the following lemma 
from [17] which studies the sensitivity of the stationary 
distribution of Markov chains: 

Lemma 2: Consider two Markov chains with transition 
matrices Q and R, and stationary distributions % and ft, re- 
spectively. We define G = I — Q, and denote its pseudoinverse 
by G 1 " = [gjj]. Hence, 

P-tf|L<M#-<2lL d2) 



where K\ 



4 



is called the condition number of the 



max,-. 7 

chain described by the transition matrix Q. 

In the next subsections, apply the above lemma to bound 
the factor ||vi (R) — vi (Q) |L. In the first subsection we com- 
pute an explicit expression for the norm of the perturbation 
\\R ~ GIL as a function of the properties of the nodes. In the 
second subsection, we study the coefficient K s in terms of 
the eigenvalues of EW*. 



A. Infinity Norm of the Perturbation \\R — GIL 

Our approach is based on studying the entries of the n 2 x 
n 2 matrix R = E [Wk <g> W^], and compare them with the entries 
of Q = EW k <$EW k . The entries of Q and R are of the form 
E(wy)E(Hy s ) and E(w, 7 w rl ), respectively, with i,j,r, and s 
ranging from 1 to n. These entries can be classified into 
seven different cases depending on the relations between the 
indices. In the Appendices II and EH, we present explicit 
expressions for each one of these cases. A key observation 
in our approach is to notice the pattern that the above entries 
induces in the matrices R and Q. For sake of clarity, we 
illustrate this pattern for n = 3 in Fig. 1, where the numbers 
in parenthesis correspond to each one of the seven cases 
identified in the Appendices. 

Since the entries of R follow the same pattern as the entries 
of Q (although these entries are different), the perturbation 



matrix A = R — Q also follows the same pattern as R and 
Q. Hence, comparing the entries of Q and R we can easily 
deduce the following seven cases for the entries of A: 



A 2 
Ai 

A 3 



A 5 = A 7 = 0, 

■M 2 (zi)-Mj ( Zi ) 
a 



^(-M 2 ( Zi )+M 2 (zi)) 



A 4 = a; 



Mi(zt)-M 2 (zi) (l-Mife-)) 2 



di 



A 6 



1+2V/-3M, 
di(di-l) 



tjia si (l-Mi(zi)y, 



d 2 



for dj> 1 , 
for d,- = 1 . 



where Mj = E [l/(J;+ l) 2 ] can be written as a hyperge- 
ometric function that depends on p, and di (see Appendix 
III). From the above entries, we can compute via simple, but 
tedious, algebraic manipulations the following expression for 
the infinity norm of the perturbation: 



|*-GL = ||A|| 



= max {Si} , 

!</<« 



(13) 



where 



2(1 



M 



M 



Note that 5, is a function that depends exclusively on the 
sequence of nodes properties {(pu,d„)} ueV ; hence, \\R — <2|L 
depends exclusively on the set of nodes degrees and probabil- 
ities, but it is independent on how these nodes interconnect. 
In the following subsection, we show how the pattern of 
interconnection among nodes influences the upper-bound of 
the variance in (TTOb via the condition number K s in (fT2l . 

B. Perturbation-Based Bound for the Variance 

In this subsection we derive an explicit relationship be- 
tween the condition number K s and the network structure 
via the spectral properties of EW^ =£+(/— £)D~ . This 
result will then be used to bound the variance of x* in ( TTOb . 
We base our analysis on a the following bound, derived by 
Meyer in [18], relating K s with the eigenvalues of Q, denoted 

as {Mi}"=i : 



max 



< 



2(n 2 



L 



(14) 



s(l-Wt) 

Before we present our main result, we need some notation. 

2 

Denote by {A,}" =1 and {l~lj} j =l the set of eigenvalues of EWt- 
and EWjigiEWfc, respectively. The ordering of the eigenvalue 
sequences is determined by their distance to 1, i.e., 1 1 — Aj| < 
1 1 — Xj\ for i < j. Hence, our result can be stated as follows: 
Theorem 3: The variance of the asymptotic consensus 
value of <[TJ can be upper-bounded by 



var(x*)<||x(0)®x(0)||. 



(A) 



max {Si} 

l<i<n 

1b) 



n(i 




(C) 



(15) 



where {A,-}" =1 are the eigenvalues of EWj, and the product 

n(i - kkj) = H(ij) s .t. (i - ■ 

Proof: We start our proof from (TTOb 

var(x*) < ||vi(/?)-v 1 (e)m||x(0)®x(0)|| 1 

<K\\R-Q\\J\x(0)®x(0)\\i 
^Ksfnmx {S,}) \\x(0)e>x(p)\\i 

\l<i<n / 

(c) 2(« 2 - 1) 

< n i „ \ m ^ H x (°) 0x Wlli ' 
IK 1 — MfcJ 

where we have used Lemma [2] in inequality (a), the expres- 
sion for — filloo m 03 m equality (b), and the upper 
bound in dT4b in step (c). We obtain the statement of the 
theorem by applying the following standard property of the 

2 

Kronecker product: {li k }" k=l = {XjXj} ^. . <;( . ■ 

Remark 4: The bound in ( TTBT l separates the variance into 
three multiplicative terms representing each one of the fol- 
lowing elements (for convenience, we have underlined each 
one of these terms in $15\): 

(A) This first term exclusively depends on the initial 
condition as indicated by ( TTTb . 

(B) The second term depends solely on the nodes 
properties c/, and pf. 

(C) The last term represents the influence from the 
overall graph structure via the spectral properties 
of EW k . 

Remark 5: It is specially interesting to study the implica- 
tions of Term (C) in the asymptotic variance. For example, 
given the sequences of degrees and probabilities, {dj}" =l 
and {pi}" = i, the influence of the network structure on the 
variance is given via Term (C). Since the eigenvalues of 
EW k are key in this term, it is interesting to briefly describe 
the homogeneous Markov chain with transition matrix P = 
EW k =£+(/ — E) D- X A T C . This Markov chain presents a self- 
loop in each one of its n states with transition probability 
Pa = M- , as well as a link from ; to j with transition 
probability P i} = fl-Mf'jdr 1 for £ E c . From an 
analysis proposed by Meyer in [18], we conclude that Term 
(C) is primarily governed by how close the subdominant 
eigenvalues of EW k are from 1. In particular, the further 
the subdominant eigenvalues of EW k are from 1, the lower 
the upper bound in dT5Y In the next section, we illustrate 
the relationship between spectral properties of EW k and the 
variance of the asymptotic consensus value with several 
numerical simulations. 

V. Numerical Simulations 

In this section we present several numerical simulations 
illustrating our results. In the first subsection, we numerically 
verify the result in Theorem Q] In the second subsection 
we present some examples to illustrate the implications of 
Theorem[3]in the variance of the asymptotic consensus value. 



j k 




Fig. 2. The above figure represents a chain of leader nodes, denoted as 1,2, 
and 3, with a set of followers. Connections between leaders are represented 
as pairs of directed links with existence probabilities p m , and p vu , with 
u,v 6 {1,2,3}. We have also included in this figure the connection of each 
leader and one of its followers (drawn as a smaller circle on top of each 
leader). In this case, the directed link incoming the follower (leader) exists 
with probability p/f (p/;). The rest of followers are represented under each 
leader, where we represent each pair of directed links as a single undirected 
edge for clarity. 

A. Expectation of the Asymptotic Consensus 

In our first simulation we take a graph G c composed by 3 
stars connected in a chain (see Fig. 2). This graph is intended 
to represent, in a social context, three leaders with a set of 
followers. In Fig. 2, we represent the leaders using large 
circles marked as 1,2, and 3, and the followers using smaller 
circles. We assume that each follower only listen to one 
leader (the center of a star) and nobody else. In this particular 
example, the first, second and third leaders have 4, 8, and 16 
followers, respectively. In each time step, a directed random 
graph is built by choosing a set of directed communication 
links from G c . We have fixed the probability of existence of a 
directed link coming into followers to be equal to p/f — 1 /2, 
for all such links (see Fig. 2). Also, the probability of 
existence of a directed link coming into a leader is inversely 
proportional to the degree of the leader. More specifically, 
we have chosen p m . = l/d v , for v= 1,2, and 3, where p m , 
represents the existence of a directed communication link 
from the «-th node to the v-th leader. 

In this example, we compute the asymptotic consensus 
value for 100 realizations of a random consensus algorithm 
with initial conditions Xi(0) = i/N. We represent the his- 
togram of these realizations in Fig. 3, where the empirical 
average of the 100 realizations equals 0.5077. The theoretical 
expectation for the asymptotic consensus value applying 
Theorem □ equals E(x*) = ZiPWi{i/N) =0.4595, which is 
in great accordance with the empirical value. For future 
reference, we have also computed the empirical standard 
deviation of the 100 realizations to be 0.0298, and the three 
eigenvalues of EW k closest to 1 are {0.9823,0.9449,0.75}. 

B. Expectation of the Asymptotic Variance 

In this subsection, we numerically illustrate some of the 
implications of Theorem [3] in the variance of the asymptotic 
consensus value. Although the upper bound stated in this 
Theorem is not tight, there are some qualitative implications 




Fig. 3. Empirical histogram for 100 realizations of the random directed 
consensus over the graph in Fig. 2. The empirical average and the empirical 
standard deviation of the realizations are 0.4595 and 0.0298, respectively. 




Fig. 4. Ring of 3 leaders and corresponding followers. 



VI. Conclusions and Future Work 



that are consistent throughout our numerical experiments. 
For example, as we mentioned in Remark [5] given a set of 
initial conditions and nodes properties, the upper bound in 
(T5[ is primarily governed by how close the subdominant 
eigenvalues of EWj are from 1. In particular, the further 
the subdominant eigenvalues of EWk are from 1, the lower 
the upper bound. We illustrate in the following numerical 
experiments that the lower the upper bound, the lower one 
should expect the variance to be. 

In the first experiment, we slightly modify the network 
described in the previous subsection and study the influence 
of this modification on the eigenvalues of EWj and the 
variance of x*. Our first modification is a change in the 
probabilities of existence of a directed link from the u- 
th node to the v-th leader without changing the network 
topology. In particular, we choose the new probabilities to be 
Puv = 3/c/y, for v = 1,2, and 3. In the modified network, the 
eigenvalues of EW k closest to 1 are {0.9789,0.9372,0.75}. 
Hence, since the effect of our modification on the eigenvalues 
is to move them away from 1, we should expect the variance 
of x* to be reduced according to Remark [3] In fact, running 
100 random consensus algorithms with the same initial 
conditions, xi (0) = i/N, using our new probabilities, we 
obtain a standard deviation 0.0286 (which is less than the 
0.0298 obtained before). 

In the second experiment, we illustrate how the larger 
the gap between the eigenvalues of EWk and 1, the smaller 
the variance of x*. In this case, apart from keeping the 
new set of probabilities p m , = 3/d v , we convert the 3-chain 
of leaders into a 3 -ring of leaders, as depicted in Fig. 
4. In this case, the three eigenvalues of EWk closest to 
1 are {0.9577,0.9212,0.75}, which are even further away 
from 1 than in the second example. Hence, as expected, 
the standard deviation of 100 random consensus algorithms 
is even smaller than in the second example, in particular 
0.0274. In conclusion, our simulations are consistent with 
Theorem [3] and with the qualitative behavior described in 
Remark [5] 



We have studied the asymptotic properties of the consen- 
sus value in distributed consensus algorithms over switching, 
directed random graphs. While different aspects of consensus 
algorithms over random switching networks, such as condi- 
tions for convergence and the speed of convergence, have 
been widely studied, a characterization of the distribution 
of the asymptotic consensus for general asymmetric random 
consensus algorithms remains an open problem. 

In this paper, we have derived closed-form expressions 
for the expectation of the asymptotic consensus value as a 
function of the set of initial conditions, {x u (0)} :l€V , and the 
set of nodes properties, {(p u ,d u )} ueV , as stated in Theorem 
[TJ We have also studied the variance of the asymptotic 
consensus value in terms of several elements that influence 
it, namely, (/) the initial conditions, (it) nodes properties, and 
(Hi) the network topology. In Theorem [3] we have derived 
an upper bound for the variance of the asymptotic consensus 
value that explicitly describes the influence of each one of 
these elements. We also provide an interpretation of the 
induce of the network topology on the variance in terms 
of the eigenvalues of the expected matrix EWk- From our 
analysis we conclude that, in most cases, the variance of 
x* is primarily governed by how close the subdominant 
eigenvalues of EWk are from 1 . We have checked the validity 
of our predictions with several numerical simulations. 



Appendix 
Computing the Entries of EWk 

We start by computing the entries of EWk- The diagonal 
entries of EWk are given by: 



Ew H 



1 



l+di 



k=0 K + 1 



1 di+l 



Pi (di + l) 



Similarly, the non-diagonal entries of EW^ result in: 



Similarly, we also have: 
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Appendix 
Entries of EW k ® EW A . 

We now compute the possible entries in g = EW^ ® EW^. 
The entries of the Kronecker matrix EWa- ® EWit, with EWa- = 
E+ (/ — E)D A r , present entries that can be classified into 
seven different cases depending on the relations between the 
indices. These are the cases, where we assume that all four 
indices i,j,r, and s are distinct: 



Qi=E 
Qi = E 
Q 3 =E 

g 4 = E 
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Appendix 
Entries of EW k <& Wt 

We now turn to the computation of the elements of E[Wt<S) 
Wit], which are of the form E(w, ; w„-). Again, we classify the 
entries into seven different cases depending on the relations 
between the subindices: 
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= I ^(l+2Mf)-3M) 1 )), far4>l, 
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= E(wyW r /) = E(wyWri) 
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